Could the Hercules satellite be a stellar stream in the Milky Way halo? 
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Abstract. We investigate the possibility that Hercules, a recently discovered Milky Way (MW) satellite, is a 
stellar stream in the process of formation. This hypothesis is motivated by Hercules' highly elongated shape 
as well as the measurement of a tentative radial velocity gradient along its body 1 1 1. The application of simple 
analytical techniques 1 2 1 on radial velocity data of its member stars provides tight constraints on the tangential 
velocity of the system (u, = — 16*| 2 km s _1 , relative to the Galactic Standard of Rest). Combining this with its large 
receding velocity (145kms~') and distance (138 kpc) yields an orbit that would have taken Hercules to within 
6+' kpc of the Galactic centre approximately 0.6 Gyr ago. This very small perigalacticon can naturally explain 
the violent tidal destruction of the dwarf galaxy in the MW's gravitational potential, inducing its transformation 
into a stellar stream. 



1 Introduction 

The Hercules stellar system is one of the so-called 'ultra- 
faint' satellites of the Milky Way. Residing at a heliocentric 
distance of 138 kpc, its highly elongated morphology (el- 
lipticity of 0.65, [3 4 5 1) begs the question of how it could 
have become so distorted whilst currently inhabiting the 
outer domains of the stellar halo. The detection of a ten- 
tative radial velocity gradient in its member stars [ jj adds 
further enigma to this My - -6.6 [4| system. In this work 
(6), we investigate whether Hercules could have undergone 
a relatively recent close encounter with the Galaxy, result- 
ing in its morphological distortion from a regular dwarf 
spheroidal galaxy into the highly elongated structure ob- 
served today. 



2 Determining the orbit for Hercules 

Ascertaining the proximity of Hercules' last perigalactic 
passage requires its orbit to be calculated. In the follow- 
ing sections, we first provide a brief description of how the 
tangential velocity of an object may be determined from 
a radial velocity gradient measured along it. This is fol- 
lowed by an outline of the model used to deduce the ve- 
locity distribution of Hercules stars and the results of our 
calculations. In this study, we use a model for the Milky 
Way composed of a Miyamoto-Nagai disk and bulge Q, 
and an NFW |8| profile with parameters constrained by 
Xueetal. 0. 
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2.1 Finding the tangential velocity 

For a given Galactic potential tp, the radial velocity gradi- 
ent along a system can be expressed in terms of its helio- 
centric distance, D, radial velocity, u T , and tangential ve- 
locity, v t , as follows [2 1 Q 



V t + (V r lff) 



D 



(1) 



where x is the angle along the stream and V r i[r is the com- 
ponent of the gravitational acceleration due to the Galactic 
potential along the line of sight to the system. This results 
in the following two solutions for the tangential velocity of 
the orbit at the location of the system: 



1 

* = 2 



du, 



(2) 



For Hercules, all variables and parameters on the right- 
hand side of this equation are either known or calculable, 
hence the direct measurement of a velocity gradient over 
the body of the satellite translates to only two possible or- 
bits, corresponding to the positive and negative solutions 
of u t . 



2.2 Model and results 

We employ a maximum likelihood algorithm on Hercules 
member stars with radial velocity determinations G] to 



1 Unless stated otherwise, velocities are quoted relative to the 
Galactic Standard of Rest and distances are given as heliocentric 
quantities. 
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Fig. 1. Marginalized relative likelihood distributions for the radial 
velocity gradient, tangential velocity and perigalacticon, whose 
best values are AvJ&x = 25 ±15 km s~'/deg, Vt - -16*5, km s _I 
and ^ pcri = 6*2 kpc, respectively. From top to bottom, the hor- 
izontal, dashed lines intersect the likelihood distributions at the 
boundaries of the 1-, 2- and 3-cr confidence intervals. 



find the likelihood distribution functions of the following 
properties of the stellar system: the mean radial velocity of 
the member stars, D r , the radial velocity gradient, dvjdx, 
and the velocity dispersion of the stars, s, around the ve- 
locity gradient. Under our assumption that Hercules is a 
stellar stream, the direction of its motion on the sky is given 
by the direction of its elongation, in other words the ma- 
jor axis of the system. The position angle of the satellite 
[4 1 is therefore used as a prior in deducing the likelihood 
distributions, and hence the best model. 

Using equation (0, the likelihood distribution function 
of the radial velocity gradient can be transformed into that 
for the tangential velocity. The distance is assumed to be 
D — 138+7 kpc, obtained by averaging the distance mea- 
surements of Aden et al. ifTUl and Sand et al. [5|. Each 
value of the tangential velocity completes the 6D phase- 
space information necessary to derive an orbit and there- 
fore translates to a single pericentric distance. This leads to 
the marginalized relative likelihood distributions shown in 
Figure Q] for the radial velocity gradient, tangential veloc- 
ity and perigalacticon, with the best values and 1-cr errors^ 
for these parameters being dujdx — 25 + 15 km s~' /cleg, 
v t - -16^2 km s 1 and/? per i = 6^ kpc, respectively. These 
constitute our best model for the Hercules system under the 
assumptions and prior stated above, and the properties of 
the corresponding orbit are shown in Figure [2] Employing 
the other possible solution for the tangential velocity from 
equation (f2]i leads to orbits that never venture close enough 
to the Galactic centre to have been affected by tides, and 
are therefore considered not relevant for the purpose of this 
study. 



3 Summary 

The technique for determining the tangential motion of a 
stream via its radial velocity gradient is a powerful tool, 



2 This corresponds to the 1-cr confidence interval of the likeli- 
hood distribution, defined such that the relative likelihood drops 
by a factor e~*~ /2 at the boundaries of a k-cr confidence-interval 
range. 
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Fig. 2. Orbit of the Hercules stellar stream using our best model 
as given in the text, integrated forward and backward from the 
present location for 2 Gyr in each direction. The panels show the 
sky projection of the orbit in Galactic coordinates (top), its helio- 
centric distance (middle) and radial velocity relative to the Galac- 
tic Standard of Rest (bottom). The solid and dotted lines respec- 
tively denote the forward and backward-integrated orbit from the 
present location of Hercules, which is indicated by an arrowhead 
in the direction of the orbit. 



whose strength in aiding orbit determinations has a partic- 
ularly timely role in the era of large sky photometric and 
spectroscopic surveys. We use this technique and a max- 
imum likelihood algorithm to find a viable orbit for the 
Hercules satellite that would have taken the stellar system 
to within 6 kpc of the Galactic centre ~ 0.6 Gyr ago. Such a 
small perigalacticon could explain a disruption of the sys- 
tem so violent that we now observe it as an unbound stellar 
stream. This simple possibility for Hercules should serve 
as a cautionary warning for treating all recently discovered 
ultra-faint dwarf galaxies as a single population of bound 
stellar systems. 
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